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Abstract —This paper is concerned with the state estimation 
probiem for genetic regniatory networks with time-varying deiays 
and reaction-diffusion terms under Dirichiet boundary condi¬ 
tions. It is assumed that the noniinear reguiation function is 
of the Hiii form. The purpose of this paper is to design a 
state observer to estimate the concentrations of mRNA and 
protein throngh avaiiabie measnrement outpnts. By introducing 
new integrai terms in a novei Lyapnnov-Krasovskii functionai 
and empioying Wirtinger-based integrai ineqnaiity, Wirtinger’s 
ineqnaiity, Green’s identity, convex combination approach, and 
reciprocaiiy convex combination approach, an asymptotic sta- 
biiity criterion of the error system is estabiished in terms of 
iinear matrix inequaiities (LMIs). The obtained stabiiity criterion 
depends on the upper bounds of the deiays and their derivatives. 
It shonid be highlight that if the set of LMIs are feasible, 
the desired observer exists and can be determined. Finally, two 
nnmericai examples are presented to illustrate the effectiveness 
of the proposed designed scheme. 

Index Terms —Genetic regulatory networks. Reaction-diffusion 
terms. State estimation, Wirtinger-based integral ineqnaiity. 


I. Introduction 

I n the last decades, due to the increasing progress in genome 
sequencing and gene recognition, genetic regulatory net¬ 
works (GRNs) have become a significant area in biological 
and biomedical sciences. However, there still exists large gap 
between the genome sequencing and the understanding of 
gene functions which have become challenge problems in 
system biology. A great amount of experimental results show 
that mathematical modeling of GRNs can be a powerful tool 
for researching the gene regulation process and discovering 
complex structure of a biological organism in-El. Gener¬ 
ally, there are two basic models for GRNs: Boolean model 
(discrete-time model) a, a and differential equation model 
(continuous-time model) 0-181. Differential equation model 
describes the change rates of the concentrations of mRNAs and 
proteins. Furthermore, differential equation model has been 
most frequently utilized since it can more precisely describe 
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the whole network and made it possible to understand the 
dynamic behavior of whole network in detail. 

In biological systems, particularly in GRNs, stability is 
the most significant and essential dynamical behaviors 0- 
Ga. It is related to not only the structure and function 
of an organism but also the strength and characteristics of 
the external disturbances. In addition, as it is well known, 
time delays caused by the slow processes of transcription 
and translation in real GRNs. It has been well shown from 
present research results that time delay may lead to instability, 
bifurcation or oscillation for systems d-nsi. However, 
mathematical modeling of GRNs without introducing delays 
will lead to wrong predictions of the concentrations of mRNAs 
and proteins. Therefore, the problem of stability analysis for 
biological systems with time delays has stirred increasing 
research interests and a great deal of excellent results has 
been reported in literature in recent decades (see, e.g., 0, 

0, d-EH). 

In some mathematical modeling, it is implicitly assumed 
that the genetic regulatory systems are spatially homogeneous, 
namely, the concentrations of mRNA and protein are homoge¬ 
nous in space at all times. However, there are some situations 
in which these assumptions are not reasonable. For instance, 
it might be necessary to consider the diffusion of regulatory 
proteins from one compartment to another 0, m-m. 
In this situation, the general functional differential equation 
model can not precisely describe genetic regulatory process 
more or less. Hence, it is imperative to introduce reaction- 
diffusion terms in mathematical modeling of GRNs. To the 
best of authors’ knowledge, the delayed GRNs with reaction- 
diffusion terms are only studied in Il25]l - ll28l . Ma et al. lIZTll 
introduced reaction-diffusion terms to GRNs for the first time 
and established delay-dependent asymptotic stability criteria. 
Based on the Lyapunov functional method, Zhou, Xu and Shen 
E&i investigated finite-time robust stochastic stability criteria 
for uncertain GRNs with time-varying delays and reaction- 
diffusion terms. Han and Zhang ||25]| . gradually improved 
Ma et al.’s results by introducing novel Lyapunov-Krasovskii 
functional and employing Jensen’s inequality, Wirtinger’s in¬ 
equality, Green’s identity, convex combination approach and 
reciprocally convex combination (RCC) approach. 

In complex biological networks such as neural networks 
and GRNs, it is often the case that only partial information 
about the states of the nodes is available in the network 
outputs. In order to understand biological networks better, it 
is indispensable to estimate the states of the nodes through 
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available measurements. Hence, the problem of state estima¬ 
tion for biological networks has been one of the investigated 
dynamical behaviors in recent years a, m-M- However, 
to the best of the authors’ knowledge, there is still no any 
published result on the state estimation problem for GRNs 
with time-varying delays and reaction-diffusion terms, which 
arouses our research interests. 

Motivated by above discussion, we aim to investigate the 
state estimation problem for GRNs with time-varying delays 
and reaction-diffusion terms. By introducing new integral 
terms into a novel Lyapunov-Krasovskii functional and em¬ 
ploying Wirtinger-based integral inequality, Wirtinger’s in¬ 
equality, Green’s identity, convex combination approach and 
RCC approach, an asymptotic stability criterion of the error 
system is established in terms of LMIs. Thereby, a state 
observer is designed, and the observer gain matrices are 
described in terms of the solution to a set of LMIs. 

The rest of the paper is organized as follows; the problem 
is formulated and some preliminaries are given in Section 2; 
in Section 3, an asymptotic stability criterion for the error 
system is established, and an approach to design state observer 
is proposed; two numerical examples are provided in Section 
4; and finally, we conclude this paper in Section 5. 

Notation We now set some standard notations, which will 
be used in the rest of the paper. I is the identity matrix with ap¬ 
propriate dimension, represents the transpose of the matrix 
A. For real symmetric matrices X and Y, X > Y{X > Y) 
means that X—Y is positive definite (positive semi-definite). H 
is a compact set in the vector space ffi." with smooth boundary 
9H. Let C^{X,Y) be the Banach space of functions which 
map X into Y and have continuous fc-order derivatives. For a 
positive integer n, let (n) be the set {1, 2 ,..., n}. 


H. Model description and preliminaries 


This paper considers the following GRNs with time-varying 
delays and reaction-diffusion terms 1251 : 


' dmi{t,x) I 


dt 

dpi{t,x) 

dt 


=ELi 


dxk 


Ytik 


dmi{t, x) 
dxk 


- a^fhiit.x) 


+ E?=i - a{t), x)) + Qi, 




d 


D 


dpi{t,x) 


ik ' 


dxk 


- CiPi{t,x) 


+birhi{t - T{t),x), 


( 1 ) 


where i £ (n), x = col(a;i, a:2,..., x/) € ft C ft = 
< Lk,k £ (()}, Lk is a given constant; D^k > 0 
and D*i^ > 0 denote the diffusion rate matrices; rhi{t,x) and 
Pi(t,x) are the concentrations of mRNA and protein of the 
ith node, respectively; Ui and Ci are degradation rates of the 
mRNA and protein, respectively; bi represents the translation 
rate; W := \wij\ £ is the coupling matrix of the genetic 

networks, which is defined as follows: 


! 7ij, if j is an activator of gene i, 

0, if there is no link from gene j to 
—liji if j is a repressor of gene i, 

here 7^^ is the dimensionless transcriptional rate of transcrip¬ 
tion factor j to gene i; gj represents the feedback regulation 


function of protein on transcription, which is the monotonic 
function in Hill form, i.e., gj{s) = where H is the Hill 

coefficient; qi = Ii is the set of all the nodes which 

are repressors of gene i; a{t) and T{t) are time-varying delays 
satisfying 

0 < r(f) < T, T(f) < /ii, 

0 < u{t) < CT, ait) < p2, 

where r, a, pi and p 2 are non-negative real numbers. 

The initial conditions associated with GRN ([U are given as 
follows; 

rhi{s,x) = (j)i{s,x), x £ s £ [—d, 0], 2 £ (n), 

Pi{s,x) = cj)*{s,x), X € H., s £ [—d, 0 ], 2 £ (n), 

where d = max{CT,T}, and (j)i{s,x), (()*{s,x) £ C'^([—d, 0] x 

H,K). 

In this paper, the following type of boundary conditions 
(Dirichlet boundary conditions) is considered; 

cc) = 0, a: £ OH, t £ [—d, -boo), 

Pi{t,x) = 0 , a: £ OH, f £ [—d,-boo). 

Now, we assume that 


m*{x) := col(m*(a:), r7i2 (x),... , m* (x)) 


and 

p*{x) := co\{p*i{x),p* 2 {x),... ,p*(x)) 

are the unique equilibrium solution of GRN ([T]i, that is. 


0= ELi 


0 

dxk 


Dk 


dxk 


— aim* (x) 


+ + 9i, 


0= ELi 


0 

dxk 


D. 


.dp*{x) 

dxk 


- Ctp* (x) + bim* (x) 


for i £ (n). Obviously, the transformations, fhi = rhi—m* and 
Pi = Pi — P*, i ^ {n)! transform GRN ([T]i into the following 
matrix form: 

Odt(f,x) _ i f dfh{t,x) \ 

dt ~ ^^==1 dxk V " Oxfc j 

-Am{t,x)+ Wf{p{t-a{t),x)), 
dp{t,x) _ I 0 f^^ dp{t,x) \ 

dt ^^==1 dXk \ dXk J 

—Cp{t,x)+Bm{t — T{t),x), 


where 


A = diag(ai, 02,..., a„), C = diag(ci, C2,..., c„), 

B = diag(&i,&2, • ■ -An), 

-Hfc lfiag(fdl/i; 7 Pd 2 ki • ■ • 7 P^nk): 
Dl=diag{D*i,,D;,,...,D:,), 
m{t, x) = col(mi(<, x), ffi2{t, x),..., mn{t, x)), 
pit,x) = CO\{pi{t,x),P2{t,x),.. . ,Pn{t,x)), 
f{p{t - a{t),x)) 

= C 0 l(/i(pi(f - cr(f),x)), • • ■ ,fn{Pu{t - a{t),x))), 
Mpiit-a{t),x)) = gi{p^it-a{t),x)+p*)-gi{p*), i £ (n). 
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Because of the complexity of GRN ([3]), it is normally of 
the case that only partial information about the states of the 
nodes is available in the network outputs. In order to obtain the 
true state of Q, it becomes necessary to estimate the states 
of the nodes through network measurements. The available 
measurements are given as follows: 


Zp(t,x) = Np(t,x), 


dp{t, x) 


+Wf{p{t - a(t),x)) 
+Ki[zm{t,x) - Mm(t,x)], 


d 


^'==1 dxk 


dp{t, x) 


- Cp{t, x) 


dt 


D* 

dxk V ^ 9xk 
-{C + K 2 N)p{t, x) + Bm{t - T{t),x), 


where 


/.(O) = 0,0 < ^ < e R,t/ 3 ^ 0,i e (n), 

y 


namely. 


(4) 


where Zm{t,x) and Zp{t,x) are the actual measurement 
outputs, and M and N are known constant matrices with 
appropriate dimensions. 

To estimate the states of GRN Q through available mea¬ 
surement outputs in (01, we construct the following state 
observer: 


We introduce the following lemmas which play key roles in 
obtaining the main results of this paper. 

Lemma 1 (Jensen’s Inequality): 041 . 051 For any con¬ 
stant matrix = M > 0 of appropriate dimension, 

any scalars a and b with a < b, and a vector function 
w : [a,b] —>■ K" such that the integrals concerned are well 
defined, then the following inequality holds: 

f J u>(s)ds j M f yu;(s)ds j < {b—a) jvF{s)Mw{s)ds, 


{la le w{s)dsd9'^ M {j^^ Jg w{s)dsd0^ 


< 


(b-a)^ rb rb T 


la Iff w'^{s)Mw{s)dsd9. 


dt - UXk \ UXk 

+Bm{t - T{t), x) + K 2 [zp(t, x) - Np{t, x)], 

(5) 

where rh{t,x) and p{t,x) are the estimations of m{t,x) and 
p{t,x), respectively, and Ki and K 2 are the observer gain 
matrices to be designed later. 

The initial conditions for the state observer (|3]) are assumed 
to be {m^{t,x),pi{t,x)) = {(j)i{s,x),(j)*{s,x)). 

Our aim is to find suitable observer gains Ki and K 2 , 
so that m(t,x) and p{t,x), respectively, approach to m(t,x) 
and p{t, x) as t ^ -foo. Let the error state vectors be 
m{t, x) = mlt, x) — m{t, x) and p{t, x) = p(t, x) — p{t, x). 
Then it follows from (01, (0i and (0 that 

'dm{t,x) _^i 9 (^ dm{t,x) 

dt a 1 k „ 

-{A + KiM)m{t, x) + Wfipit - a{t), x)), 
dp{t,x) _ I d (^^dp{t,x) 

— 2^k=l 


Lemma 2 (Wirtinger-based Integral Inequality): II36II For 
given a symmetric positive definite matrix Q G and 

a differentiable function uj : [a, b] —M", the following 
inequality holds: 


vF {u)Qw{u)du > 


1 


5- 


Oq 


Q 


LIq 

Oi 


where Q = diag(Q, SQ), LIq = w{b) — w{a) and 


= w{b) + iu(a) — 


b — a , 


w{u)du. 


Lemma 3 (Wirtinger’s Inequality): M37II Assume that the 
function / G C'^([a,6],R.") satisfies /(a) = f(b) = 0. Then 

ljHv)dv<^^^J\nv)rdv. 


Lemma 4: 11251 Let iVi > 0 and N 2 > 0 be a pair of 

diagonal matrices. Then the states of ® satisfy 


( 6 ) 


r dmJ' {s,x) ^ _d_ (p. am(t,x) \ , 

Jn dt l^k=l dxk dxk ) 




1 dxu 


7 -) d ( dm{t,x) ^ 


f{p{t - a{t), x)) = f{p{t - a{t), x)) - f{p{t - a{t),x)). 

From the relationship among fi, fi and gi, one can easily 
obtain that 


/t 






Q dt 


n* d ( dp(t,x) 

dxk I dt 


dx, 


dx. 


/(O) = 0, F izKfiz) - Kz) < 0, Vz e K", (7) 


where K = diag(5i, ^ 2 , • ■ ■, C«) > 0. 

In this paper, we assume that etTor system satisfies 
Dirichlet boundary conditions: 

mi{t,x)= 0 , X G dLl, f € [—d,-foo), 
Pi{t,x)= 0 , X G dLl, f € [—d,-foo). 


Lemma 5 (RCC Lemma): 1381 Let fi, f 2 , ■■■, Jn : D ^ . 
have positive finite values, where D is open subset of 
Then the RCC of fi over D satisfies 

niin|jj^.Q,^>Q ai=i} X/i 
= Ei /* (t) + maxg,_^ (t) gi,j (t) 


subject to 

g,, : R™ 


— 9 i,j ( 0 ; 


hit) 9iFt) 
9t,jit) hit) 


> 0 . 
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III. Observer design 

In this section, we will design a state observer © for GRN 
that is, find a pair of observer gain matrices Ki and K 2 
such that the trivial solution of system I© is asymptotically 
stable under Dirichlet boundary conditions. For this end, we 
define 

~ 0l4nxn5 


Ci COl(0j.^ (i—l)n7 Onx ^ 


ip{t,s,x) = col ( m(s,x) 


m(s, a;)ds 


tf}{t,s,x) = col ^{s,x), J p{s,x)dsj, 

<,{t,x) =co\{m{t,x),m{t — T,x),m{t — T{t),x),p{t,x) 
p{t - a, x),p{t - a{t), x), f{p{t, x)), 
f-(p(t-ait),x)),^^,^, 

/t-r(t) 2;)ds, 

^ //_,(*) p{s, a^)ds, = 3 ^ p{s, x)ds). 

Theorem 1: For given scalars r, ct, pi and p 2 satisfying 
0, the trivial solution of error system ® under Dirichlet 
boundary conditions is asymptotically stable if there exist 
matrices Qf = Qi > 0 (i € (5)), Rj = Rj > 0 {j £ (4)), 
MJ = Mj > 0 {j £ (2)), diagonal matrices Pj > 0, Aj > 0 
(j £ (2)), and matrices Gi, G 2 , Wi and W 2 of appropriate 
sizes, such that the following LMIs hold for r £ {0, f} and 
cr £ {0, if}: 


Rj 


R, G, 

GJ R, j 


> 0 , 


j G (2), 


( 8 ) 


^{T,a) = $o + $i + $2(T,CT) + $3 + $4(T,CT) + $5(r,CT) < 0, 

(9) 

where 

$0 = — 2e'iAiej + e^^AiK + e^KAie'^ — 2esA2e^ 

+egiFA 2 e|' + eg,A2Ke^ — 69 (^ 1 ^ + WiM)e^ 
—ei{PiA + WiM)^eg +e9PiVFe|' 

+e8VF^Pie^ - 269^16^ - eio(P2G + W2N)eJ 
-ei{P2G + W2Nfelo + eioP2Pe|’ 

+e3P'^P2e}o - 2eioP2efo, 


$3 = -O.STT^eiPiPie}’ - 2ei(PiA + WiM)e'^ 
PeiPilFeg PesW^PieJ 
-0ATT‘^eiP2DleJ - 2ei{P2G + W2N)eJ 
+e4P2Pe|’ + e3P^P2ef, 


^’ 2 (r,CT)= eiQief - (1 -/ri)e 3 Qie|’ 

+ 64(3364 - (1 - /i 2 ) 66(3366 

+ Ai( 32 Ai' + t[AiQ 2 A^ + A 2 Q 2 A.f) 

— A 3 Q 2 A^ — t{A 2 ,Q 2 A 2 + A 2 Q 2 A'^) 
PA^Q 2 A'^ + AqQ 2 A'^ 

+r(A 5 Q 2 A|’ + A6Q2AI’) + 0 iQ 40 }’ 

+(t( 0 iQ 40 I’ + 02Q40D - 
—criQ^Q^Qr^ + 02 ( 3403 ^) + 

PQ^QaQI + <j[Q^Qa &1 + 06 Q 40 D. 

d ’3 = e-iQb^j — (1 — ^2)68(3568, 

$ 4 (t, cr)= $41 - $ 42 (r) - $43(0-) 

— [Ay A8]Pi[A7 A 8 ]'^ — [07 08 ]P 2[07 08 ]”^, 

O '“P O P* O 'P O 

^41 — ^ egi^iCg + G eioi? 26 io “ 1 “ ^ O' 641^46^ , 

$ 42 (t) = f(f - r)ei 2 P 36 f 2 + freiiR^eii, 

$43(t) = d{d - a)euRie{^ + d-tTei3P4ef3, 
$5(GCr)= $51 — $52 — $53 

-(I^AsMiAi^ - e^08M20i’, 

^2 

$51 = —egMiel + —ewM 2 elQ, 

$ 52 = (61 - eii)Mi(ei - eii)^ 

+ (e3 — 6i2)Mi(e3 — 612)^, 

$53 = (64 — 6i3)M2(64 — 613)^ 

+ (66 — 6 i 4 )M 2(66 — 644)’^, 

Ai = [61 'f6i2], A2 = [60 611 — 612], 

A3 = [62 fei2], A4 = [fei2 f^6i2], 

A5 = [611 — 612 r(6ii — 612)], Ag = [60 61 — 62], 

Ay = [63 — 62 63 + 62 — 2ei2], 

Ag = [61 — 63 Cl + 63 — 2eii], 

01 = [64 (fei 4 ], 02 = [60 613 — 644], 

03 = [65 0-644], 04 = [d-644 0-^644], 

05 = [643 — 644 d-(643 — 644)], 0 g = [60 64 — 65], 

07 = [6g — 65 6g + 65 — 2644], 

08 = [64 — 6g 64 + 6g — 2643], 

P4 = diag(P 4 , 3 P 4 ), R2 = diag(P 2 , 3 P 2 ), 

Ml = 4diag(M4, 3 M 4 ), M 2 = idiag(M2, 3 M 2 ), 

T a 
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D*l = diag 




I 





i 


-.E 

fc=i 


Ll) 


and Lfe, Di^, D*f,, A, B, C, W and K are the same with 
previous ones. 

Moreover, the observer gain matrices are given by Ki = 
P^^Wi and K 2 = Pi^W 2 . 

Proof: Construct a Lyapunov-Krasovskii functional for 
error system (|6]l as follows: 


5 

V{t,m,p) = 

2 = 1 


where 


Vi{t,m,p) 


= {t,x)Pim{t,x)dx 

+ jQP^(t,x)P2p{t,x)dx 

I 7-, dm{t,x) 

Jn 

+eU/. 


dxk 

dp^{t,x) 


n dxi. 


-0-dx 

OXk 

* ap(£^ 1 

dXk 


V 2 {t,m,p) = mA{s,x)Qim{s,x)dsdx 

+ fa /E x)Q 2 p(t, s, a:)dsda; 

+ In ft-a(t) ^)Q3 P{s, a::)dsdx 

+ fn ft- - s, x)Qii'ip(t, s, a:)dsdx, 


V 3 {t,m,p)= / (j}{s,x))Q 5 f{p{s,x))dsdx, 

JQ J t — 




dm-^{s,x) R^djn^dsdOdx 


t -\-6 ds ^ ds 

+lL 

+ T /o !-r ft+e m'^{s, x)R3m{s, a;)dsd0d2: 

In f-w ft +0 a:)i? 4 p(s, a;)dsd0da:, 


Vnit, m,p)=J„ S"_, f° /.VK^m, 2dte2d,dAdfldx 


Taking the time derivatives of Vi{t,m,p) {i G (5)) along the 
trajectory of error system (|6]l yields 


■§^V2{t,m,p) 

= J^m'^{t,x)Qim{t,x)dx 

— (1 — f(t)) mf'it — T(t),x)Qim{t — T{t), x)dx 

+ jQp'^{t,x)Q3p{t,x)dx 

-(1 - a{t)) J^p'^it - a{t),x)Q 3 p{t - a{t),x)dx 

+ J^p'^{t,t,x)Q2p{t,t,x)dx 

— -T,x)Q 2 p{t,t -T,x)dx 

+2 fn ft-T x)Q 2 ^^^%f^dsdx 

+ fn t, x)Q 4 'tp{t, t, x)dx 

— t — a, x)Q 4 'ilj{t, t — a, x)dx, 

+2 fn fLw x)Q 4 ^'^^1'^’'"^ dsdx 

^ fn ^)‘^ 2 (T(t),cr(t))(;(t, x)dx, 

( 11 ) 

§jV3{t,m,p) 


= -(1 - F{p{t - a{t), x))Q 3 f{p{t - a{t),x))dx 

+ fn x))Q5f{p{t, x))dx 


dt 


\W,,n.p)= T=/„a!!2!dafl,2=!faldi 


( 12 ) 




.2 f 


dt 


r rt drfd^R M£A1 Asdx 
"JnJt-a ds ds 

+r^ /n x)R3m{t, x)dx 

-rfnt _ mf'{s, x)R3m{s, tcjdsdx 
+CT^ Fp^{t, x)R4p{t, x)dx 
~^fn ftf^p'^(s,x)R4p(s,x)dsdx, 


dt 


v,{t. „.,p) = i /„ ad^asiMja^dx 


(13) 


- /n S^dpdSdp 

75-2 r 

+~/r 


Jt +9 

dp^ (t,x) , j dp(t,x) 
fn- 2 A+ 


da: 




ds 


Msd^da;. 


(14) 

From Green formula, Dirichlet boundary conditions and 
Lemma 13 we have 


-§^Vi{t,m,p) 

— {A + KiM)m{t, x) + Wf{p{t — a{t),x)) 

+2 4p^(i,x)pJEL,!fe(d>:.^) 


d24„pm,x)p,[Y;,,,^AD;s%A) 

— {C + K 2 N)p{t, x) + Bm{t — T(t), x) dx 

+2EL./„ a^^p.Dt3|-(a:fdi)d* 

lOV^* r dp'^(t,x) □ n* d t dp(t,x) ^^^ 
'X^2^k=lJn dxk ^^^kdxk^ dt 


dx 


( 10 ) 


= 2 EL./»("i^((.»)r.i)»a 5 g£i)'_, 

-2Ei../„a^p.D.d^dx 

= -2EL./„a:S^P.D^d^dx 


• n ds 


< 


^z^fc=iJn dxk dxk 

jQrn^{t,x)PiDL'm{t,x)dx, 
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where 


(m'^{t,x)PiDk 


dm{t,x) 

dxk 

dm{t,x) 




Similarly, 

2ELi/nP^(^>3;)P2^ 

^ InP'^it^^)P2Dlp{t,x)dx. 

The combination of (fTOl l. (fTSl l and (fTSl l gives 

— {A + KiM)m{t, x) + WJ{p(t — (T(f), x))] dx 
+2/aP^(^:a;)P2 -^Dlp{t,x) 

— {C + K 2 N)p{t, x) + Bm{t — T{t),x)] dx 

+2ELi In 

r dp'^{t,A r> n* d i dp(t,x) -^^^ 

2 ^k=l Jn dxk dxk •' dt 

Q idm(t,x)^ 


On the other hand, by Lemma [T] it yields 

—T rri^{s,x)R3m{s,x)dsdx 

= —T TO^(s, x)R3m{s, x)dsdx 




<-W^ In ILt'^ a:)dsP3 to(s, x)dsdx 
-Tmfn ft-T(t) a;)dsP3 x)dsdx 


t-r(t) _ 


(16) 


r(t) Jn Jt-T(t) 

^ “ /n x)$ 42 (T(f))<t(f, x)dx. 

In the same way, 

-^/n/t-7fP^('S’2;)-R4P(s ,x)dsdx 
^ -/a?’^(i,a^)®43(o-(f))<7(f,x)dx. 

Combining (fOl l. (l20l l- (l2Tt yields 


( 22 ) 


(23) 


dm (t,x) p p g ^ am(t,x) -^^ 

r dp^{t,x) p n* 3 t dp(t,x) 

' 2-^k—l Jn fixu 3xu ^ 


+^T,LJn^-^P^D 


= r Jn dxk ^‘2^kdxk'^ ' dt ')^^- 


(17) 

Note that the second term on the right of (fTSl) can be written 


(25) 


as; 


c ft dm^{s,x) p dm{s,x) n„. 

^JnJt-T ds ds risux 

— -T ^ r«-2~(*) dm'^(s,x) P dm{s,x) ,, 

f ft dm'^{s,x) p dm(s,x) ,, 

^JnJt-r(t) ds ds 

Applying Lemma |2] one can obtain that 


(18) 


< 


( ft—■-(«) dm'^(s,x) p dm(s,x) n 1 

^ Jn Jt-T ds _ 1 ds nsux 

x) A7P1 Af <j(<, x)dx 




and 


f ft dm^(s,x) n 9 m(£^ , , 

^JnJt-T(t) as ds risux 

^ J^)^8-RrA8^(f,x)dx. 

This, together with ([8]l, (fTST l and Lemma |5] implies that 

f ft dm^(s,x) p am(s,x) . , 

' Jn Jt-T 9s 9s risux 

^ “/n^^(^’^)[^7 A8 ]Pi[A 7 A8]^<r(f,x)dx. 

Similarly, 

- “/ n ^^(^’^)[®7 08 ] P 2[07 08 ]^^( C 2 ;) dx . 


(19) 


( 20 ) 


■§fV4{t,m,p) < ^^^(f,x)$4(T(f),cr(f))<?(f,x)dx. 

(24) 

The second term on the right of (fT4l i can be divided into 
three parts: 

2=3:EM,a^d.d«d* 

= A. ^^“■^'I.dfldx 

- In r;“’ /«"’ ^^=£^d.d«d. 

-(T - -(()) /„ »=!Ei£)M,»=Md,dx. 

By using Lemma [1] we can estimate the following inequalities 

- In I°r,„ 2=SH4M,»=^d.d«dx 

- In /.‘+7“’ 2=4?^M, »;t^d,dWx 

< -Thin Hi., /.'« S=4ll^id,dW,TAt, 

x/-„.)C^'l,d9dx 

-wAwin /:;*'* 2=4?^d,d«dxAt, 

x/„/-V'‘*/.+7'”^“«dx 

= -2/j^<j'^(J,x)$52<?(f,x)dx. 

(26) 

As in (fTOl l. the last term on the right of (l25l l can be bounded 
by applying the same procedure, 

-(7 - r(l)) /„ 

< (j'^(f,x)A8MiA|’^(f,x)dx. 

(27) 

In a similar manner, 

- In I°,m A. ^2S^Aft^d,dMx 

- /□ /.+7'‘' a!^Aft2£feEd»d9dx (28) 


- -2/n7^(i,2;)4'537(f,x)dx, 


and 


-(n- T(t)) a!^Af,Sl!^d,dx 

< - J^(^J’(i^a;)08M20|’<?(f,x)dx. 


( 21 ) 


(29) 



























































7 


Combining (fT^ and (|25])-(|29l), we can obtain 
d f 

-^V5{t,m,p) < j ^^(i,x)$5(r(i),cr(t))<j(t,a:)dx. (30) 

It is also easy to see that 


2 r dm^(t,x) p 

Jq dt t 


Sfe=i offe ( 


D 


dm{t,x) 


(A + KiM)m{t, x) + Wf{p{t — a{t),x) 


(31) 


d7n{t,x) 

dt 


dx = 0 


and 


2/. 


ap'^(t,x 
n dt 


■Po 


9 f T^* Op{t,x) ^ 
2^k=l dxk \ k dxk ) 


(C + K 2 N)p{t, x) + Bm{t — T{t),x) 
dx = 0. 


(32) 


dp{t,x) 

dt 


According to Lemma|4] Green formula and Dirichlet boundary 
conditions, we have 


r) f dm^(t,x 

^ Jn 9t 

.T! 


■P, EL. 2k (P-^) 


= ^Inrn^it 2 x)PiEk=i^ 


D 


d / dm{t,x) 


= -2EU/t 

Similarly, 


dm^ {t,x) 


=iJn^kJ-PiDk^r-W^)dx. 


k dxk dt 

9 ! dm(t,x) ^ 
dt 


dx (33) 


9 r dp^(t,x) p 9 / 

^ Jq dt ^ 2 -^k—l dxk y 


D 


« dp(t,x) 


da; 


= 2/„p^(i,i)PjELi»fe 


k dxk 

n* 9 / dp(t,x) 

^k dxk y 9t 


dx 


(34) 

Finally, for the diagonal matrices Ai > 0 and A 2 > 0, it 
can be obtained from (|7]i that 

2f'^{p{t,x))Ai_f{p{t,x)) 

- 2p'^{t,x)KAif{p{t,x)) <0, 


(35) 

(36) 


2/^(p(f - o-(f),x))A 2 /(p(< - a{t),x)) 

— 2 p^{t — a{t), x)KA 2 f{p{t — a{t), x)) < 0. 

From (fTTT l. (fT2l l. (fTTl l, (l24l l and (l30ll-(l36ll, one can obtain 


< /n‘^^(^i3;)$(T(f),cr(f))<r(f,x)dx. 

Since ^{T{t),a{t)) depends affinely on rit) and ait), respec¬ 
tively, it follows from (|9|l that ■^V{t,m,p) < 0 for all T{t) 
and a{t) satisfying (|2]l. Therefore, the trivial solution of error 
system (|6]l is asymptotically stable. This completes the proof. 


We end the section by the following remarks on Theorem 


ffl 

Remark 1: Compared with 
paper are as follows: 

1) We introduce new integral items like 


the advantages of this 


n j —T Jt-\-9 


m'^{s, x)R^m{s, x)dsd0dx 


and 


0 /.O 


* dm'^(s,x) dm{s,x) , 

^ 'Ml -^2_J.dsdAd6»dx 


T J 9 J t+A 


ds 


ds 


into Lyapunov-Krasovskii functional and employ Wirt- 
inger-based integral inequality (instead of Jensen’s in¬ 
equality) to estimate the derivative of the second one, 
which will get more accurate result. 

2) The so-called convex combination approach and RCC 
approach are employed simultaneously, which will im¬ 
prove the precision of estimation to the concentrations 
of mRNA and protein. 

3) The coefficients of some items in <j(f,x), like and 

, play a very important role in simplification of 
the LMI condition (|9|l. 

4) We use ip{t,s,x) and 'ip{t,s,x) instead of m{s,x) and 
p{s,x) in V 2 {t,m,p), respectively. This will highly 
maintain consistent with V^{t,m,p). 

Remark 2: The approach proposed in this paper can easily 
be applied to establish a delay-dependent and delay-rate- 
dependent asymptotic stability criterion for GRN ([T]l. Due to 
Remark [T] above, the criterion is certainly less conservative 
than ones in Ha, EH, El. 


IV. Illustrative examples 
I n this section, two numerical examples are provided to 
demonstrate the effectiveness and applicability of the proposed 
state observer. 

Example 1: Consider GRN El* with measurements El, the 
deterministic parameters are given as: 


A = diag(0.2,1.1,1.2), B = diag(1.0,0.4,0.7), 


C = diag(0.3,0.7,1.3), Li = L 2 = L 3 = 1, 


W = 


0 

-0.5 

0 


0 

0 

-0.5 


-0.5 

0 

0 


Di = D 2 = D 3 = diag(0.1,0.1,0.1), 
D* = D* = D* = diag(0.2,0.2,0.2), 


M = 

0.5 

-0.6 

0 

0.3 

0.8 

-0.2 

N = 

0.7 

-0.25 

0.3 

0.4 

0.2 

-0.3 


Here the regulation function is taken as fix) = One 

can get © holds when K = 0.657. When t = a = 3 and 
pi = P 2 = 2, by using the MATLAB YALMIP Toolbox, one 
can see that the LMIs given in Theorem[T]are feasible with the 
following feasible solution matrices. To save space, we only 
list some of the feasible solution matrices as follows: 


Pi = diag(57.6506,44.1104,50.5774), 


P 2 = diag(25.7909,39.4682,32.9357), 




1.3165 

-0.0077 

-0.0077 

Qi — 


-0.0077 

2.2460 

0.0907 



-0.0077 

0.0907 

1.5577 



■ 2.8401 

0.1585 

0.0103 ■ 

Qb 


0.1585 

5.6978 

0.0748 



0.0103 

0.0748 

3.3839 








































0.3 > 


Ri = 

R 2 = 

Ml = 

Wi 


5.8721 

0.0099 

0.0606 

2.0575 

-0.0035 

0.0052 

1.6422 

0.0046 

0.0301 


0.0099 

3.7715 

-0.1517 

-0.0035 

3.6837 

- 0.1101 

0.0046 

1.4934 

-0.0324 


0.0606 

-0.1517 

4.2376 

0.0052 

- 0.1101 

2.8021 

0.0301 

-0.0324 

1.2873 


34.7528 

9.8144 

4.3837 


25.0882 

-15.2841 

9.3021 



Figure 1. The real trajectory of mRNA (m(t,x)) 


W2 = 


16.3285 18.6193 

5.3137 -9.3968 

-12.5043 22.4536 


Moreover, we can get the corresponding observer gain matri¬ 
ces as follows; 


Ki = Pf^Wi = 
K 2 = P2^W2 = 


0.6028 

0.4352 

0.2225 

-0.3465 

0.0867 

0.1839 

0.6331 

0.7219 

0.1346 

-0.2381 

-0.3797 

0.6817 



Example 2: When I = n = 1, GRN (O is simplified into 


dm{t^x) 

dt 


dp{t,x) _ 

dt 


dx 


Di^^)-Arn{t,x) 


Figure 2. The estimated trajectory of mRNA (rh(f, x)) 


+Wf{p(t-a{t),x)), 


+Bfh{t — T{t),x). 


_a_ 

dx 


(37) 


We choose the values of parameters in (ITtI i are as follows, 


A = 0.2, B = l, C = 0.3, Li = 1, 

W = -0.5, Di = 0.1, Dl = 0.2, 

M = 0, N = 0.7. 

When pi = = ‘2, K = 0.65 and t = a = 1, 

for Dirichlet boundary conditions, by using the MATLAB 
YALMIP Toolbox to solve the LMIs given in Theorem [U we Figure 3. The estimation error of mRNA (m(4, x) - rh{t, x)) 
obtain the following feasible solution matrices. To save space, 
we only list some of the feasible matrices as follows; 



Pi = 1.8102, Qi = 0.0120, 


i?i = 0.2847, Ml = 0.1556, 

Wi =0, W2 = 1.1478. 

Moreover, we can get the corresponding observer gain matri¬ 
ces as follows; 

Ki = Pf^Wi = 0 , 

K2 = P2^W2 = 1.5571. 

Further, when a{t) = T{t) = 1, the state responses of GRN 
Ozli, observer © and the corresponding error system are given 
in Figures [T}© 



Figure 4. The real trajectory of protein (pit, x)) 
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= 0-4^ 


I 0 . 3 - 



0 0 


Figure 5. The estimated trajectory of protein (p{t, x)) 


0.3^ 



0 0 


Figure 6. The estimation error of protein (p(t, x) — pit, x)) 


V. Conclusions 

In this paper, the state estimation problem for a class of 
GRNs with time-varying delays and reaction-diffusion terms 
are studied. An state observer is designed to estimate the gene 
states through available sensor measurements, and guarantee 
that the error system is asymptotically stable. By introducing 
new integral terms in a novel Lyapunov-Krasovskii functional 
and employing the so-called Wirtinger-based integral inequal¬ 
ity, Wirtinger’s inequality. Green’s second identity and, convex 
combination approach, RCC approach, a sufficient condition 
guaranteeing the existence of state observers is established 
in terms of LMIs. The concrete expression of the desired 
state observer has been presented in Theorem [T] Finally, two 
numerical examples are given to illustrate the effectiveness of 
the theoretical results. 
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